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Abstract — This paper considers the distribution of the mu- 
tual information of frequency-selective spatially-uncorrelated 
Rayleigh fading MIMO channels. Results are presented for 
OFDM-based spatial multiplexing. New exact closed-form expres- 
sions are derived for the variance of the mutual information. 
In contrast to previous results, our new expressions apply for 
systems with both arbitrary numbers of antennas and arbitrary- 
length channels. Simplified expressions are also presented for 
high and low SNR regimes. The analytical variance results 
are used to provide accurate analytical approximations for the 
distribution of the mutual information and the outage capacity. 

Index Terms — MIMO Systems, Orthogonal Frequency Divi- 
sion Multiplexing, Mutual Information 



I. Introduction 

Multiple-input multiple-output (MIMO) antenna technology 
has emerged as an effective technique for significantly improv- 
ing the capacity of wireless communication systems. A great 
deal of work has been done on analyzing the MIMO capacity 
in various flat-fading channel scenarios, since the pioneering 
work of [1] and [2]. In particular, the mean (ergodic) capacity 
has now been comprehensively investigated (e.g. see [3-18] 
and references therein). 

In addition, the outage capacity has also been investigated 
for flat-fading channels. This is an important capacity mea- 
sure for systems with stringent delay constraints, and also 
provides information about the system diversity [19]. With 
the exception of the exact two/three antenna results presented 
in [20,21], outage capacity analysis has typically involved 
approximating the distribution of the mutual information, 
since exact closed-form solutions are not forthcoming. It has 
been shown that the Gaussian distribution provides a good 
approximation in many cases [5,8, 13,22,23]. 
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TABLE I 

MIMO-OFDM PRESENCE IN STANDARDS 



Standard 


Technology 


WLAN IEEE 802. 1 In 


OFDM 


WiMAX IEEE 802.16-2004 


OFDM/OFDMA 


WiMAX IEEE 802. 16e 


OFDMA 


MBWA IEEE 802.20 


OFDM 


WRAN IEEE 802.22 


OFDM 


3GPP Release 8 


OFDMA 



In this paper, we consider frequency-selective MIMO chan- 
nels, which are applicable for many current high data-rate 
wireless systems. We focus on MIMO orthogonal frequency- 
division multiplexing (OFDM) systems, since they form the 
underlying technology for a many emerging MIMO standards, 
as summarized in Table I, and consider spatial multiplexing 
transmission. Despite their key practical significance however, 
for these systems (and indeed frequency-selective MIMO 
channels in general) there are relatively few analytic MIMO 
capacity results. The ergodic capacity (average mutual infor- 
mation) was considered in [19,24—26] and [27,28], assuming 
Rayleigh and Rician channels respectively, and was found 
to be easily obtained by summing the equivalent flat-fading 
ergodic MIMO capacity of each individual OFDM subcarrier 
In contrast, the outage capacity does not decompose in this 
way. 

Calculating the outage capacity for frequency-selective 
channels is difficult due to the non-negligible correlations 
between subcarrier channel matrices. As such, the investi- 
gation of outage capacity has usually been performed using 
simulation studies [19,29,30]. It appears that the only cur- 
rent analytical outage capacity results for frequency-selective 
channels are presented in [31], [32] and [33], all of which 
derive a Gaussian approximation for the mutual information 
distribution. The results in [31] however, are based on deriving 
exact expressions for the mutual information variance of 
single-input single-output (SISO) channels only; whereas the 
results in [32] and [33] are based on approximating the mutual 
information variance using asymptotic methods. Specifically, 
[32] considers multiple-input single-output (MISO) channels 
with asymptotically large channel lengths, whereas [33] con- 
siders MIMO channels with infinite numbers of transmit and 
receive antennas. We note also that for the extreme frequency- 
selective fading case, ie. where the MIMO subcarrier matrices 



2 



IEEE TRANSACTIONS ON INFORMATION THEORY (TO APPEAR) 



are independent across frequency, the variance of the mutual 
information could be easily calculated by adapting known 
MIMO flat-fading variance results given, for example, in [8] 
and [15]. For many practical systems however, the subcarrier 
channels are typically highly correlated across frequency, and 
this approach cannot be applied. 

In this paper, we consider MIMO OFDM-based spatial 
multiplexing systems with finite numbers of antennas, and 
operating over spatially-uncorrelated Rayleigh fading channels 
with finite delay spreads. We first derive new exact closed- 
form expressions for the mutual information variance. We 
also give explicit reduced formulas for the specific cases 
of multiple-input single-output (MISO), single-input multiple- 
output (SIMO), and single-input single-output (SISO) systems. 
Moreover, simplified closed-form expressions are derived for 
the variance in the high and low signal-to-noise ratio (SNR) 
regimes. 

Based on the new analytic variance results (along with 
known analytic mean results), we then present new approx- 
imations to the mutual information distribution of OFDM- 
based spatial multiplexing systems. In particular, we present a 
new closed-form Gaussian approximation, which is shown to 
be extremely accurate for many different system and channel 
scenarios. In the low SNR regime, we also present a new 
analytic Gamma approximation, which we show to be more 
accurate than the Gaussian approximation in this case. 

Finally, we use the analytic Gaussian approximation to 
estimate the outage capacity. We find that the approximation 
is very accurate, and show that for outage levels of practical 
interest, the outage capacity depends heavily on the delay 
spread of the channel. 

The paper is organized as follows. In Section II we describe 
the frequency-selective MIMO channel model, the OFDM- 
based spatial multiplexing signal model, and the associated 
mutual information. In Section III, we present the main analyt- 
ical contributions of the paper, namely, analytical expressions 
for the variance of the mutual information. The proofs are 
relegated to the appendices. In Section IV, we approximate 
the distribution of the mutual information, and investigate the 
outage capacity. 

The following notation is used throughout this paper. Ma- 
trices are represented with uppercase boldface, and vectors 
with lowercase boldface. The superscripts {■)'^ ,{■)*, and {-Y 
indicate matrix transpose, complex conjugate, and complex 
conjugate transpose respectively. The matrix Ip denotes apxp 
identity matrix. We use dct (•) and tr(-) to represent the matrix 
determinant and trace operations respectively. The operator 
E[-] denotes expectation, and Var( ) denotes variance. The 
real Gaussian distribution with zero-mean and unit-variance 
is denoted Af{0, 1), the corresponding complex circularly 
symmetric Gaussian distribution is denoted CAf{0, 1), and the 
chi-square distribution with r degrees of freedom is denoted 

Xr* 

II. OFDM-Based Spatial Multiplexing Systems 

A. Channel and Signal Model 

We consider a single-user OFDM-based spatial multiplexing 
system employing Nt transmit antennas, Nj. receive antennas. 



and N subcarriers. The channel is assumed to be frequency- 
selective and is modeled as a length-L finite impulse-response 
(FIR) filter (as in [19,25]), for which the discrete-time input- 
output relation is given by [25] 

L-l 

y[q] = ^pH[p]x[q -p]+ n[q] (1) 
p=0 

where yi[q\ £ C^*^^ is the signal vector transmitted at sample 
index q, y[q] £ qN^-xi jg fj^g corresponding received signal 
vector, and ■n.\q\ £ C^'^^ is the noise vector containing inde- 
pendent elements ^ 0^/(0, 1). Also, ap, for p = 0, . . . , L — 1, 
represents the channel power delay profile, and is normalized 
according to 

L-l 

E^' = i- (2) 

p=0 

The Nr X Nt random matrices H[p], for p = 0, . . . , L — 1, 
represent the MIMO channel impulse response. These matrices 
are assumed to be mutually uncorrected, and are assumed to 
be known perfectly at the receiver but are unknown at the 
transmitter.. The channel is assumed to be quasi-static, re- 
maining constant for the duration of a codeword, but changing 
independently from codeword to codeword. Throughout the 
paper, we assume that the channel elements exhibit spatially- 
uncorrelated Rayleigh fading', in which case each H[p] con- 
tains independent elements ^ C7V(0, 1). 

At the transmitter, the time-domain input sequence x[q] is 
generated as Nt parallel OFDM symbols. The symbols for 
each antenna are OFDM modulated using an A^-point inverse 
fast-Fourier transform (IFFT) prior to transmission. At the 
receiver, OFDM demodulation is performed at each receive 
antenna using an A^ -point FFT. A key advantage of OFDM- 
based spatial multiplexing is that equalization is simple, since 
the frequency-selective MIMO channel is transformed into A^ 
orthogonal flat-fading MIMO subchannels via the IFFT/FFT 
processing. 

To maintain orthogonality in the presence of intersymbol 
interference caused by multipath, OFDM systems typically 
employ a cyclic prefix extension. Assuming that the cyclic 
prefix is longer than the delay spread of the channel, we can 
write the equivalent frequency domain input-output model for 
OFDM-based spatial multiplexing as follows 

rfe = Hfcafe + nfe, fc = 0,...,A^-l (3) 

where a^^ is the transmitted vector for the fcth subcarrier, as- 
sumed to be i.i.d. Gaussian with covariance matrix E [afea[,] = 
^^Iatj, Tfc is the received vector for the fcth subcarrier, and 
Uk is the corresponding complex AWGN vector satisfying 
i?[nfen|i] — lN^6[k — i], where S[-] is the Kronecker-delta 
function. Also, H^, is the fcth subcarrier channel matrix given 



Note that a number of recent investigations liave studied the impact of 
spatial correlation on MIMO capacity (see, eg. [8,9,12,13]). We do not 
follow this line of work here however, since our primary focus is to study the 
impact of frequency-selective fading on capacity, in which case the effect of 
correlation is observed across frequency. 
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by 



L-l 

Hfc = ^ fjp U[p\ exp ( -~j2TT—p 

p=0 



(4) 



containing independent entries {'H.k)ij ~ CJ\f{0, 1). Note that 
due to the finite-length impulse response, correlation exists 
between different subcarrier channel matrices. Using (4), the 
correlation coefficients between the channel elements on two 
arbitrary subcarriers k and £ is easily derived as follows (see 
also [34]) 

-'^"/""dit-mj-f] (5) 



Pk-e 



E 

p=0 



^2 ^-j2.ik-e 



for all As expected, these frequency correlation 

coefficients depend only on the difference between subcarriers 
(i.e. k — £), and not on the subcarriers themselves. 

Note that with the above model, the SNR per receive 
antenna per subcarrier (henceforth referred to as 'the SNR') 
is given by 7. 

B. Mutual Information 

The focus of this paper is on the statistics of the mutual 
information of OFDM-based spatial multiplexing systems. It 
is now well-known that the instantaneous mutual information 
in b/s/Hz for a given channel realization is given by [19] 



N ^' 

k=0 



(6) 



where Tk is the instantaneous mutual information for the fcth 
OFDM subcarrier, given by 



log2 det liv. 



J_ 
Nt 



(7) 



Note that the loss in mutual information due to the cyclic 
prefix has been neglected in (7). The mean (ergodic) mutual 
information is given by 

1 

E [Jotdm] = ^ E ^ [^fe] ■ 
fc=0 

It is obvious that (8) is equivalent to the ergodic mutual 
information of a flat-faded channel, for which case closed- 
form expressions are now available [6, 15,35]. 

III. Variance of the Mutual Information 

In this section we derive new closed-form expressions 
for the variance of the mutual information of OFDM-based 
spatial multiplexing. Our results are exact, and apply for 
arbitrary finite system and channel parameters. We also present 
simplified expressions for the variance in the high and low 
SNR regimes, and give explicit reduced variance expressions 
for the cases of MISO, SIMO, and SISO systems. These results 
will be subsequently used in Section IV for providing accurate 
approximations to the mutual information distribution, and to 
the outage capacity. 



A. Exact Analysis at All SNRs 

The following theorem presents an exact expression for the 
variance of the mutual information of MIMO-OFDM systems. 

Theorem 1: The variance of the mutual information of 
MIMO-OFDM systems is given by 



det(B,,,) (Er=idet(A,))' 



TV 



Vra{n)Vm{rn) 



(9) 



where m = mm{Nr,Nt),n = nia.x{Nr , Nt) , T,n{') is the 
complex multivariate gamma function defined as 



and 



{n)^Y[T{n-i+l) 

i=l 



(E" 1 dot(A, 



Er=iEr=ie^^'/^det(a,.(p,)) 



(10) 



\Pd\ = 
< \pd\ < 1 

\pd\ = 1 
(11) 



(A.) 



2J 



(12) 



The matrix A,, is to x m, with j)'^ element 

bl for j ^ r 

ble^^/^gi{b+l) for j = r 

The matrices Br,s and Cr,s(-) are mxm with (i, j)*^ elements 
given by (13) and (14) respectively (at the top of the next 
page). Also, b — n + m — i—j, t ^ n^m, z = r + i + j — I, 
u = T + i + t, V = T + j + t, and G3 4(-) is the Meijer-G 
function (see [36, eq. (9.301)] for definition). 



5i(z) = ^E, 



h=l 



7 



92 



h=l 



Nt 



7(l-|PdP) 



(15) 



(16) 



where E/j(-) is the Exponential Integral (see [37, eq. (5.1.12)] 
for definition). The function riij{-, •) is defined as 

- 1\ / 1 - 
t 



mAfi^):pd)^nT + j)J2 



\pd\ 



X (T+j~t),_l/(z-t) 



(17) 



for an arbitrary input function /, and (•),. is the Pochammer 
symbol 



(a)r = a ■ (a + 1) 



. _ T{a + r) . . 

r(a) 



1 . 

(18) 

Proof: See Appendix I. □ 
Note that the exact variance expression in Theorem 1 can be 
easily evaluated since it primarily involves simple polynomial 
and exponential terms, as well as standard functions such as 
exponential integrals and Meijer-G functions, both of which 
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bl 

6!e^'/'^gi(5+ 1) 
2(iV,/7)''+^e^'/^Eto (?)(-!) 

6,t-6 



xg: 



,T ; |t-&.t 

^t/7|o, ■ 



t-f)-l.t-b-l,t-fc-l 



, for j ^ r and j =/= s 

, for J = r or j — s, and r ^ s 

for , j ~ r ~ s 



(13) 



(Cr,s(/3d)) 



(l-lpdl ) g^d-lp^l 



) |pd|"r(M)r(t-)g2(«)g2(^) 

=0 t!(T+i)! 



for i 7^ r, j s 

for i — r, j ^ s 

for i ^ r, j — s 

for i = r, j = s 



(14) 



2Wt 



r(n)e 

2g^./7 (i^)" (^^*/7io":;X 



-t-i,-t-i 



for Iprfl = 
for < \pd\ < 1 
for Iprfl = 1 



(21) 



are implemented as built-in procedures in various mathemat- 
ical software packages such as Maple and Mathematica. We 
also note that although Theorem 1 involves infinite series of 
exponential integrals, its numerical evaluation can be made 
more efficient by exploiting the following recurrence relations 
[37, eqs. (5.1.7) and (5.1.14)] 



E,{z) = -Ei(-z) 
E„+i{z)= - {e-'- - zEn{z)) 



(19) 



for 2: > 0. As such, only a single exponential integral must 
be explicitly evaluated when summing these series. Moreover, 
it turns out that this infinite series converges quickly, and can 
generally be evaluated with less than 20 terms. Therefore the 
computational challenge associated with this series is very low. 

The following corollary presents an exact variance expres- 
sion for the mutual information of SIMO and MISO OFDM 
systems (i.e. cases with m = l,n > 1). To the best of our 
knowledge, this result is also new. 

Corollary 1: The variance of the mutual information of 
SIMO/MISO-OFDM systems is given by 



Var(Xofdm) = 



(log2(e))^ 



T{n) 



2 

d=l 

N 



7 



E 



n-l 
t 



:;4 {Nth\ 



0,-t-l,-t-l.-t-l 



(20) 

where (p{pd) is defined in (21) at the top of the page. 

□ 

The following corollary presents an exact variance expres- 
sion for the mutual information of SISO OFDM systems (i.e. 
m — l,n — 1). 

Corollary 2: The variance of the mutual information of 



SISO-OFDM systems is given by 



Var(Xotdm) = (log2(e))^ 
2ei/'' 



JV-l 



7V2 



Y,{N^dMp,) 



d=l 



where 



viPd) 



o2/7 



(5i(l))^ 



for \pd\ = 



(1 - |p^|2) e^o^u^ 



E^olP<il'*(ff2(l + i))^ forO<|p,|<l 



%^Gti (1/7I0 for \Pd\ 



= 1 



(23) 

□ 

Very recently, an equivalent expression for the SISO-OFDM 
variance given in (22) was presented in [31]^. In contrast to 
(22) however, the equivalent result from [31] is not expressed 
in closed-form, and it requires the evaluation of infinite series 
of incomplete gamma functions. 

In Fig. 1 we compare the analytical variance expression 
(9) with the variance obtained via Monte-Carlo simulation. 
Results are presented for two different Nt x Nr antenna 
configurations as a function of the channel length L. A 
uniform power delay profile is assumed (i.e. cr^ = for 
p = 0, . . . L — 1), 8 subcarriers are used (simply by way of 
example; similar results are obtained for higher numbers of 
subcarriers), and the SNR is set to lOdB. In all cases we see 
a precise agreement between the simulated and analytic curves. 
Moreover, the variance is seen to be largest for the system with 
the least antennas, regardless of the channel length. For both 
antenna configurations, we see that the variance reduces with 
increasing L, and that this reduction is most significant for 

^Note that this expression was not exphcitly stated in [31]. It can however 
be trivially obtained by following the derivation of (47) and using [31, Eqs. 
(12), (41), and (48)]. 
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- ln(l - \pd\^)fj,A^,pd) + \pd\^(^'~^HvjAH{z - -ln(l - \pa\')vjA^, Pd))) (27) 




Channel Length (L) Channel Length (L) 



Fig. 1. Variance of the mutual information of MIMO-OFDM for different 
Nt X Nr antenna configurations, and different cliannel lengtlis (uniform power 
delay profile). The "Analytic Variance" curves are based on (9). 8 subcarriers 
is considered, with SNR of 10 dB. 



Fig. 2. Variance of the mutual information of MIMO-OFDM for different 
SNRs, and different channel lengths (uniform power delay profile). The 
"Analytic Variance" curves are based on (9). A 2 X 2 system is considered 
with 8 subcarriers. 



small L. For example, by increasing the channel length from 
L = 1 (flat-fading) to L = 2, the variance for both antenna 
configurations is more than halved. 

In Fig. 2 we plot the analytical variance expression (9) 
and Monte-Carlo simulation results for different SNRs, as a 
function of L. Again we see a precise agreement between 
the analytical and simulated results. From this figure we see 
that for a given channel length the variance of the mutual 
information varies monotonically with the SNR. This increase 
is most significant for small values of L. 



B. Analysis at High SNR 

The following theorem presents a closed-form expression 
for the variance of the mutual information of MIMO-OFDM 
in the high SNR regime. This result is simpler than the exact 
general variance result given in Theorem 1, as it does not 
involve any infinite series. 

Theorem 2: In the high SNR regime, the variance of the 
mutual information of MIMO-OFDM systems is given by 

/ 2 



^ m— 1 



N 



1 /^^^^ 



t=o 



ip{n-t) 



where 



Apd) = 



Er=iEr=idct(c.,3(Prf)) 

r™(n)r„(in) 

TAA,'A{n-t) 

+ (E'/Io^ VX^i - t)) for \pd 



for \pd\ = 
for < \pd\ < 1 



1 



(25) 



where Cr,s{Pd) is an to x m matrix with (i, j)"^ element for 
the cases i ^ r or j ^ s given by 



(Cr.s{pd)j 



Vi,j{^,Pd) for ij^r,j^s 

Vi,j {ip{z),Pd) for i = r, j s 

IPdP^'"^^ Vj,r {Az),Pd) for i^r, j = s 



(26) 



and for the case {i — r, j = s) by (27) at the top of the 
page, where i' — max(i,j) and j' = min(z,j). Also, rji,j{-) 
is defined in (17) in Theorem 1, ^. (•) is defined in (126), h{-) 
is given by 

MPd)-ln(l-|pdn-/C, (28) 

and JC — 0.5772 ... is the Euler-Mascheroni constant. The 
function H{-) denotes the harmonic number 



(24) 



H{z) 







for z > 
for z — Q 



(29) 
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n-l H{b-1) 



V{Pd) 



-1 1 

= 1 6 



ln(l-|p.n-E 



b 

b-2 
t=0 



b-t-1 



for \pd\ = 

for < \pd\ < 1 
for \pd\ = 1 



(32) 



and ip{-) is the digamma function defined as [37, eq. (6.3.2)] 



ip{n-t) = H{n-t- 1) -/C 



(30) 



with first derivative corresponding to the poly gamma 

function [37, eq. (6.4.1)]. 

Proof: See Appendix III. □ 
The following two corollaries present very simple high SNR 

variance expressions for the special case of SIMO/MISO and 

SISO systems respectively. 

Corollary 3: The variance of the mutual information of 

SIMO/MISO-OFDM systems at high SNR is given by 



Var°°(Jofdm) 
= (log2(e))' 




d) V'iPd) 



^P'{n) 
N 



(31) 



where (fipd) is given by (32) at the top of the page, with Li^ (•) 
denoting the dilogarithm function [37, eq. (27.7.1)] . □ 
Corollary 4: The variance of the mutual information of 
SISO-OFDM systems at high SNR is given by 



Var°°(Iofdm) 
= (log2(e))' 



/ o ^-1 2 \ 



(33) 

□ 

It is important to note that the results in Theorem 2 and 
Corollaries 3 and 4 do not depend on the SNR. Therefore, a 
main insight which we can draw from these expressions is 
that the variance of the MIMO-OFDM mutual information 
converges to a deterministic limit as the SNR increases, 
which we have now quantified precisely. This phenomenon is 
illustrated in Fig. 3, where we plot the variance of the MIMO- 
OFDM mutual information for different Nt x Nr antenna 
configurations, and for different SNRs. The "Analytic Variance 
(High SNR)" dashed lines are based on (24) for the 2 x 3 
case, (31) for the 1x2 case, and (33) for the 1x1 case. The 
"Analytic Variance (Exact)" curves are based on (9) for the 
2x3 case, (20) for the 1x2 case, and (22) for the 1x1 case. 
Monte-Carlo simulated variance curves are also presented for 
further verification. We see that the results converge quickly 
in all cases. 



C. Analysis at Low SNR 

The following theorem presents a very simple closed-form 
expression for the variance of the mutual information of 
MIMO-OFDM in the low SNR regime. 



Monte-Carlo Variance 

■ Analytio Variance (High SNR) 

■ Analytio Variance (Exact) 




Fig. 3. Variance of tlie mutual information of MIMO-OFDM for different 
Nt x Nr antenna configurations, and for different SNRs. The "Analytic 
Variance (High SNR)" lines are based on (24) for the 2x3 case, (31) for 
the 1x2 case, and (33) for the 1x1 case. The "Analytic Variance (Exact)" 
curves are based on (9) for the 2x3 case, (20) for the 1x2 case, and (22) 
for the 1x1 case. 16 subcarriers are considered, and the channel follows an 
8-path uniform power delay profile. 



Theorem 3: In the low SNR regime, the variance of the 
mutual information of MIMO-OFDM systems is given by 

Var°(Xofd„.) = ilog,{e)rl-f- 1 + 2 ^ 



N -d. 



\Pd\ 



NNt \ ^ N 

\ d=l / 

(34) 

Proof: See Appendix IV. □ 
The following corollary gives upper and lower bounds (as 
a function of the frequency correlation coefficients) for the 
variance of the MIMO-OFDM mutual information in the low 
SNR regime. 

Corollary 5: In the low SNR regime, the variance of the 
mutual information of MIMO-OFDM systems satisfies 

^ ^ Var°(Tofdm) ' ' 
N - Var°(Xfiat) 

where Var"(Xfiat ) denotes the mutual information variance for 
an i.i.d. flat-fading Rayleigh MIMO channel. The left-hand 
side is an equality for \pd\ = (independent fading across all 
frequency subcarriers), and the right-hand side is an equality 
for \pd\ — 1 (identical fading across all subcarriers, i.e. flat- 
fading). 
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Proof: The proof follows by using 



< 



\Pd\ 



< 1 



in (34), and noting that 



Var°(Xflat) = (log2(e))2^, 



(36) 

(37) 

□ 



which is found by directly setting = 1 in (34). 

It is interesting to note from (35) that in the low SNR 
regime, the scaling of the MIMO-OFDM variance with respect 
to the flat-fading variance depends only on the channel delay 
profile, and is independent of the number of transmit and 
receive antennas. 

For the particular case of a uniform power delay profile 
(i.e. with dp = 1/L for all p = 0, . . . , L — 1), we can obtain a 
simple insightful expression for the variance ratio in (35), as 
given below. 

Corollary 6: For a uniform power delay profile, (35) be- 
comes 

Var°(Iofdm) 



1 

N 



— < 



Var°(Xfiat) 




< 1 
(38) 

where the left-hand side is an equality for L = N, and the 
right-hand side is an equality for the case L — 1. 

Proof: The proof follows trivially from (35) after noting that 
the frequency correlation-coefficients (5) in this case can be 
expressed as [38] 



sm 



Pd 



V N ) 



Lsin(f) 



eJTT . (39) 



□ 

The summation in (38) is of a similar type to that in [33, eq. 
(60)], which gave an asymptotic expression for the variance for 
large antenna numbers, and involved the same squared-ratio 
terms. As mentioned in [33], as L increases, the ratio becomes 
more peaked as a function of d, thereby decreasing the overall 
sum. Thus, from (38) we see that the variance of the mutual 
information varies inversely with the channel delay spread in 
the low SNR regime. This agrees with previous observations 
seen via simulation studies in [19], and for the regime of large 
antenna numbers in [33]. These results are further corroborated 
in Fig. 4, where (38) is plotted as a function of the channel 
length L. 



IV. Outage Approximation of MIMO-OFDM Based 
Spatial-Multiplexing 

We now use the analytic expressions from the previous 
section to present and investigate approximations for the distri- 
bution of mutual information. We then use the approximations 
to estimate outage capacity. 

Unless otherwise stated, for all results in this section we 
model the channel according to the exponential power delay 



Monie-Carlo Variance Ratio 
■ Analytic Variance Ratio 



Flat-fading Variance Ratio 



" independent-fading Varia 



Channel Length (L) 



Fig. 4. Ratio of the MIMO-OFDM mutual information variance and the flat- 
fading MIMO mutual information variance at low SNR, for different channel 
lengths (uniform power delay profile). The "Analytic Variance Ratio" curve is 
based on (38). A 2 x 2 system at —25 dB is considered with 32 subcarriers. 
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Mutual Information (b/s/Hz) 

Fig. 5. Pd.f. of MIMO-OFDM mutual information for different A^t X TV,, 
antenna configurations. The "Analytic Gaus. Approx." curves are based on 
the exact mean formula (48) and exact variance formula (9). 64 subcamers is 
considered, with SNR of 20 dB. The channel follows an 8-path exponential 
power delay profile with Xexp = 4. 



profile [39] 







for < p < i 
otherwise 



(40) 



where iiTcxp is a parameter which characterizes the rate of 
decay of the power delay profile as a function of p, and is 
loosely related to the rms delay spread [39]. 



A. Gaussian and Gamma Approximations 

We first investigate the accuracy of a Gaussian approxima- 
tion for various system configurations and channel scenarios. 

Fig. 5 presents the analytical Gaussian approximation for the 
MIMO-OFDM mutual information p.d.f. based on the exact 
mean and variance expressions in (48) and (1) respectively, as 
well as empirically generated p.d.f.s (Monte-Carlo histogram). 
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Fig. 6. P.d.f. of MIMO-OFDM mutual information for different rms delay 
spreads (channels follow an 8-path exponential power delay profile, with 
different i^cxp)- The "Analytic Gaus. Approx." curves are based on the 
exact mean formula (48) and exact variance formula (9). 16 subcarriers is 
considered, with SNR of 20 dB. 



Fig. 7. P.d.f. of MIMO-OFDM mutual information. The "Analytic Gaus. 
Approx." curve is based on the exact mean formula (48) and exact variance 
formula (9). The "Asymptotic Gaus. Approx. (From [33])" curve is based on 
[33, eqs. (59) and (60)]. 2x2 antennas and 32 subcarriers are considered, 
with 20 dB SNR. The channel follows an 8-path uniform power delay profile. 



for different antenna configurations. A 64-subcarrier system is 
considered with SNR of 20 dB. We see that the analytic curves 
match the true distribution almost perfectly for both antenna 
configurations. We also present curves for a simulation based 
Gaussian approximation (based on the mean and variance of 
the IVIonte-Carlo generated histograms) for further verification. 
Note that these curves are indistinguishable from our new 
analytical Gaussian approximation curves. 

Fig. 6 compares the analytical Gaussian approximation with 
empirically-generated p.d.f. curves, for different channel rms 
delay spreads. Again we see that the analytic Gaussian approx- 
imation is accurate in all cases. Moreover, we see a significant 
reduction in the variance of the mutual information as the rms 
delay spread increases (i.e. as K^xp increases). Again note that 
the Monte-Carlo Gaussian approximation is indistinguishable 
from our new analytical Gaussian approximation curves. 

Fig. 7 compares our new analytic Gaussian approxima- 
tion with the asymptotic Gaussian approximation previously 
derived in [33]; formally derived under the assumption of 
asymptotically large antenna numbers. To our knowledge, this 
is the only other comparable analytical result in the literature 
which applies for arbitrary-length frequency-selective MIMO 
channels. In the figure, we consider a 2 x 2 system at 20 
dB SNR. The channel has a uniform power-delay profile; for 
which simple approximations to the mean and variance of 
the mutual information were explicitly presented in [33, eqs. 
(59) and (60)]. Clearly, although the approximation in [33] 
was shown to be quite accurate for some practical scenarios. 
Fig. 7 shows that our analytic Gaussian approximation is more 
accurate (although it is only shown for 32 subcarriers, the same 
observation has been made for all systems investigated). 

Fig. 8 presents the distribution of the mutual information at 
high SNRs, comparing MIMO, SIMO, and SISO systems. The 
analytic Gaussian approximation curves are based on a high 
SNR mean formula from [3, Theorem 2], and the high SNR 
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Fig. 8. P.d.f of MIMO-OFDM mutual information at high SNR. The 
"Analytic Gaus. Approx (High SNR)" curves are based on the high SNR 
mean formula [3, Theorem 2], and the variance formula (24) for the 2x2 
case, (31) for the 1x2 case, and (33) for the 1x1 case. 16 subcarriers is 
considered, with SNR of 35 dB. The channel follows an 8-path exponential 
power delay profile with Kcxp = 4. 



variance formula (24) for the MIMO case, (31) for the SIMO 
case, and (33) for the SISO case. We see that the analytic 
Gaussian approximation is accurate in all cases. Again note 
that the Monte-Carlo Gaussian approximation is indistinguish- 
able from our new analytical Gaussian approximation curves. 

Fig. 9 presents the distribution of the mutual information 
at low SNRs. The analytic Gaussian approximation curve is 
generated based on the low SNR mean formula obtained by 
combining (131) and (133), and the low SNR variance formula 
(34). In this case we see that a Gaussian distribution no longer 
accurately predicts the mutual information p.d.f. This can be 
explained by examining (131), where we see that at low SNRs 
the mutual information for each subcarrier is a function of 
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Fig. 9. P.d.f. of MIMO-OFDM mutual information at low SNR. The "Analytic 
Gaus. Approx (Low SNR)" and "Analytic Gamma Approx (Low SNR)" curves 
are based on the low SNR mean formula obtained by combining (131) and 
(133), and the low SNR variance formula (34). A 2 X 2 system with 32 
subcarriers is considered, with SNR of —25 dB. The channel follows a 4- 




Mutual Information (b/s/Hz) 



Fig. 10. C.d.f. of MIMO-OFDM mutual information for different rms delay 
spreads (channels follow an 8-path exponential power delay profile, with 
different Kcx_p)- The "Analytic Gaus. Approx" curves are based on the exact 
mean formula (48) and exact variance formula (9). A 2 X 2 system with 16 
subcarriers is considered, with SNR of 20 dB. 



path exponential power delay profile with 



of the channel realizations^ ie. 



which for i.i.d. Rayleigh fading is ~ X2NrNt- 
Hence, the overall mutual information (8) is distributed as 
the sum of N correlated X2NrNt random variables which (for 
small N), is clearly quite different to Gaussian. 

Motivated by this observation, we propose to approximate 
the mutual information p.d.f. at low SNR with a Gamma 
distribution. Note that a Gamma approximation was previously 
considered in the context of flat-fading channels in [40]. The 
Gamma p.d.f. is given by 

Qr^r-l -Ox 

m = % ' , X > (41) 

where r is the shape parameters and 9 is the scale parameter 
By matching the first two moments, a Gamma approximation 
for the mutual information p.d.f of MIMO-OFDM is obtained 
by evaluating 

E [lotdm] 



and 



Var (Xotdm) 



E'^ [2'ofdm] 



(42) 



(43) 



Var (lofdm) 

This analytic Gamma approximation is plotted in Fig. 9, based 
on the same low SNR analytic mean and variance formulas 
as used for the low SNR Gaussian approximation above. We 
clearly see that the Gamma approximation is much more 
accurate than the Gaussian approximation in this low SNR 
regime, and follows the simulated p.d.f. very closely. 



B. Outage Capacity 

The outage capacity lout.q is defined as the maximum 
information rate guaranteed to be supported for 100(1 — q)% 



(44) 



where q denotes the outage probability, and is thus directly 
obtained by inverting the c.d.f of Xofdm- If the distribution of 
the mutual information is Gaussian, then the outage capacity 
can be computed from the derived mean and variance as [32, 
eq. (26)] 

Iout,q = E [Xofdra] - \/Var(2'ofdm)<9^^(g) (45) 

where Q{-) is the Gaussian (J-function. 

Fig. 10 plots the outage probability for channels with dif- 
ferent rms delay spreads. The "Analytic Gaus Approx" curves 
are generated by approximating the c.d.f. in (44) as a Gaussian 
distribution, and using the exact mean and variance formulas 
in (48) and (1) respectively. Clearly this analytic Gaussian 
approximation matches closely with the empirically generated 
c.d.f. (Monte-Carlo histogram) in all cases. Moreover, we see 
that for outage probabilities of practical interest (e.g. q = 
1%), increasing the rms delay spread can yield a significant 
improvement in outage capacity. 

V. Conclusions 

This paper has considered the mutual information distri- 
bution of frequency-selective MIMO channels, in the context 
of OFDM-based spatial multiplexing systems. Exact closed- 
form expressions were presented for the mutual information 
variance, applying for arbitrary finite system and channel 
parameters. These results were used to provide accurate analyt- 
ical approximations for the distribution of mutual information, 
and the outage capacity. We observed that for most scenarios 

^Strictly speaking, computing the outage capacity would require performing 
a numerical optimization over all possible input distributions, as discussed in 
[1]. Here however, we adopt a common slight abuse of terminology, and use 
the term outage capacity to denote the outage rate for the case of OFDM-based 
spatial multiplexing systems with equal power Gaussian inputs. 
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a Gaussian approximation is accurate, while also noting that 
for low SNR a Gamma approximation yielded even higher 
accuracy. 

Appendix I 
Proof of Theorem 1 

Proof: By definition, the variance of the mutual information 
is given by 

Var(Xofdm) = E [X^^^^J - [j^^^^] . (45) 
Noting that E [lofdm] = E [Iflat], and using (6), we have 



we have 



Var(Iofdm) = E 



N-l N-1 



N2 



E^ [Xfla 



1 
1 



k=0 e=Q,i^k 



E^ [2^fiat] 



fc=0 



^ fc=o e=oj^k J 

(47) 

where Xfiat denotes the mutual information of a flat-fading 
channel. Note that the last line followed by noting that, under 
the assumptions in Section II-A, the channel statistics for each 
subcarrier (and therefore, the mutual information statistics) are 
identical [19], and moreover, these statistics are equal to that 
of a flat-fading i.i.d. Rayleigh channel. The first and second 
moments of the mutual information for flat-fading channels 
has been previously derived in terms of incomplete gamma 
functions in [15, eqs. (29) and (31)]. Using [37, eq. (6.5.9)], 
we perform some basic manipulations to express these results 
in alternative simplified forms as follows 



E [Xflat] = 



E [iLt] 



log2(e) 



(log2(e))^ 



^det(Ar) 



(48) 



m m 



^ ^ ^ r— 1 s—1 



det (B,,,) (49) 



where A^ and Br,s are defined in (12) and (13) respectively. 

The challenge is to evaluate the cross-correlation of the mu- 
tual information across frequency subcarriers E [Tkli] which, 
using (7), is given by 



E \Ik1e\ — E 



logadet [In,. + 
X log2 det I I N,. 



E 



i—l ^ 



log2 1 



J_ 



where A = {K}T=i uj_ = {wj^ 



(50) 

are the non-zero 



eigenvalues of HfcHj, and H^H^ respectively. Defining 



7 



a{x) = log2 1 



(51) 



E {TkTt_ 



E 



m m 



i=i i=i 



= EE^["(^*)"(^^)]- ^52) 
i=i j=i 

Now, to evaluate the expectations in (52), we first simplify the 
problem by exploiting the symmetry with respect to the A^s 
and u!jS. To this end, let A and oj be randomly (uniformly) 
chosen eigenvalues from A and uj respectively. Then clearly 



Pr ( A = Ai , a; = cjj ) = 



1 



(53) 



for any given i e {1, . . . , m}, j & {1, 
also write 



m}. Hence, we can 

Xi,UJ LOj) 



£;[a(A)aH] = EEPr(^ 

X E [a{X)a{uj)\X — Xi,uj = lUj] 

..mm 

= ;;^EE^["(^«)"('^j-)] ■ (54) 

i=i j=i 

where the second line follows from (53). Therefore by directly 
comparing (54) with (52) it follows that 

E [Jkli] = m^E [a{X)a{uj)] . (55) 

We point out that the simplification from (52) to (55) is par- 
ticularly important, since in order to evaluate the expectation 
in (55), clearly we only require the distribution of a pair of 
arbitrarily-selected eigenvalues, A and uj. This turns out to 
be much more convenient than dealing with the distributions 
of the individual pairs of ordered eigenvalues, ie. Xi and uij, 
required to directly evaluate (52). 



The joint p.d.f of A and uj is presented in Lemma 1 in 
Appendix II. From this lemma we see that /(A, w), and 
correspondingly E [Tkle] in (55), only depends on k and 
£ through their absolute difference, i.e. since /(A, w) only 
depends on k and £ via |pfc-^|, and from (5) 

\Pk^e\ = \Pi-k\ = \Pe~k\ ■ (56) 
Therefore the left-hand summation in (47) can be written as 



N-l N-l 

E E 

k=o e=o,e=^k 



E \Ik1t_ 



N-l 

E' 

d=l 



2 > ■XN~d)E[IoId] 



(57) 



Note that for subcarrier spacings d for which the frequency 
matrices are independent (i.e. — 0) or completely correlated 
(i.e. pd = 1), the expectations in (57) are evaluated trivially 
as 



(58) 



E [lold] - E^ [Xfiat] , 

E[I„I,]^E[ll,], 

For the case < \pd\ < I such a direct evaluation is not 
possible, and we use (55) in Lemma 1 and (67) to evaluate 
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II 



'xIt (y^!|^\/A^) a(A)a(w)dAdw 



for i — r, j 



(i-|PdP)- 



;EU(V)(i^)*^W^ for.^.,,^. 



r(r+j)|pd| 

(l-lpdl 



(60) 



r(iqf^\/A^) for i = r,j = s 

r(r+j) IPdT -A \r+i-l Y^J-Wj-n IPdl^A \* 1 for — t- 1 =/ « 

^*=o I i j {i-\p^\^ )^ (r+t)! Z _ r, J ^ s 

^ ' ' ' _L_ for z^r, j = s 



(i-iPdP) 

r(r+j)|pd| 

(l-lpdp^- 



(68) 



the expectations in (57) as follows and integrating term by term using (61), to obtain 



oo 



rri rn 

X 

r— 1 s— 1 

-m(n— 1) 



X 



2^2^det(D^,, (A,w))dAdL^ ^ i!(r + t)! 

Substituting (62), (63), and (65) into (60), we perform some 



l^"*! ^^^^det(D .) basic algebraic manipulations to write (59) as follows 



r™(n)r^(m) (1 - |/3d|2) ^^^^^^ 



(] ( \\2 2Nt/-y ™ ™ 

for < < 1, where D,, ^ is an to x to matrix with (i, j)"^ ^ * 

element defined in (60) at the top of the page. for < \Pd\ < 1- The proof is completed by substituting (66) 

and (58) into (57), and then substituting (57), (49) and (48) 
into (47) and simplifying. 

Using the identity [11] □ 

In (1 + aA) X^-'e-'^dX = T{q)e'''^b-^ ^ E;, APPENDIX II 

^ h=l ^ ^ TniNT P n F OF ArRTTR ARTT Y SfT 



(61) 



Joint P.d.f. of Arbitrarily Selected Eigenvalues 
of subcarrier matrices 



we can evaluate b(i, j) and c(i, j) in closed-form as , 7 t . \ j u u-. -i 1 . j 

^ '•'^ \ ^-J/ Lemma 1: Let A and uo be arbitrarily selected non-zero 

log2(e)e^'/T(r + j) eigenvalues of the subcarrier channel matrices H^Hj. and 

~ (\ _ |p^|2~)-j H^hJ respectively. Then the joint p.d.f. of A and uj is given 



X 

t=o 

and 



c(«,i) = 



V rn{n)V ^(w^m? (1 - l/Cdp)" 

rn m 

log2(e)e^'/^r(T + z)|pdr X ^^det(D,^,(A,^)) (67) 



r— 1 s— 1 



X 



i-1 i\ / I |2 \ * rf ) where d=k — T = n~m, and Dr,s ('^,'^) is an m x to 
( ' J ( ^^"'^ I - — (63) matrix with (z, j)*'' element given by (68) at the top of the 

, where /r( ) is the modified Bessel function of the first 

respectively. We evaluate the remaining integral, a{i,j), by '^i^id [37, eq. (9.6.10)]. 
using the power series expansion 

It{x) = y-j fc! (r + fc)! Proof: From (5), we see that HfeHj^ and H^hJ are 

fe=o ■ (frequency) correlated Wishart matrices. In [41], the joint 
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a{X,Lu, ai, Pj) for i — I, j — 1 

/- N I b{X,a„f3j) = a{X,uJJ,a^,PJ)duJj for i = 1, 1 

(Da,/3(A,W) ^.A.oo (76) 

d{ai,l3j) ^ a(Ai, Wj, a^, /3j)dAidcjj for i 7^ 1, j 7^ 1 



ordered eigenvalue density for matrices of this general form x det [ (A,tJj)5/^ ( ^T^V^V^ ) ) dAa ■■■duj„ 

was evaluated for cases where the correlation coefficient was \ \1 ^ |Prf| 

real. Extending this result to complex correlation coefficients, IPdl ™*'"^'' 

and to unordered eigenvalues, we obtain the joint eigenvalue ~ r,„(n)rm(m)(m!)2 (1 — Ip^P)™ 

density ^ y^(_]^-)pcr(a)+pcr(ffl 



Ju (X, LO) 



a 13 

m!2r„(n)r„M(l-|p,pr X / •■■ / det(a(A„u;„a.,/3,))dA....dc.„. (73) 

1 - \Pd\ J where 



X det [{X,u:,)ilr ( 1 !'|^_^|2 y^)) (69) a(A,,o.„a„/3,) = a| 



where A„j(-) is a Vandermonde determinant, defined as x e i-i^dP [ - — ^^'^^ -v/ XiUj ) . (74) 

ra 

A (A) = TT (A. — A ) = det ( A-*^) (70) Expanding the determinants, integrating term by term, and re- 

V * / forming determinants, we obtain 



Note that the extension from ordered to unordered eigenvalues 



/(A,c.) 



-7n{n— 1) 



simply involved the addition of the leading 1/m!^ factor in \ ^ > r,„(?T,)rm (m)(m!)^ (1 — IptiP)"* 
(69), whereas the extension from real to complex correlation ^ y^T^f 1)p°''(")+p'^'('^) det (A lo)\ (75) 

coefficients is trivial, and the proof is omitted. ^ ^ \ "''^ ' / 



a fj 



To evaluate (67) we mai-ginalize (69) as follows where D^^^g (•) is an m x to matrix with (i, j)th element given 

by (76) at the top of the page. Reordering rows and columns 



/(A,c.) 



yields 



L" L L'" L ■ ■ ■ det (d„,, (A, = (-i)P-(.)+Po.(« det (D„,.,, (A, c)) 

(71) (77) 



where we have let Ai = A and uti = us. We evaluate these where 
integrals by first expanding the Vandermonde determinants in 
(69) according to 



a{X,uj,iJ) for i ^ ai, j = I3i 
b{X,i,j) for i^ai, j ^ Pi 



A.(A)A„M^EMr"^"^nAr- \t:;)'^ z\tZrA 

™ (78) 
X ^(-l)P"('3) -Q (72) Applying (77) in (75) we can further simplify as follows 

where the sums are over all permutations « = {«!,..., a^} - p / (jy,\(m\\i (\ _ In.p-)'" 
and /3 = {/?!,..., of {1,...,to}, and (-I)p-W and ^ '"^ ^^^^ ^ 

(_l)pcr(/3) denote the signs of the permutations. Substituting ^ 2^ det (U„i,/3i (A, u)) 

(72) and (69) into (71) yields 

/(A,c.) 



a /3 

^2 



A2 



|p^|-m(n-l)((^_l),)^ 
r„,(n)r„(TO)(m!)2(l-|prf|2)^ 



I „ I — m(?T, — 1) 

iPdl ^ ^ mm 

r™(n)r^(m)(TO!)2(l-|p,|2)™ X ^ ^ det(D„„^,(A,w)) 



rm(n)r„(TO)TO2 (1 - |prf|2) 



nAr-^E(-i)'^"'^'^n' 

*=1 /3 j = l 
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xX^X]det(D,,,(A,c^)) . (79) 



The result now follows by combining (74), (78) and (79), 
and by evaluating the integrals c(-) and d{-) inside the 
remaining determinant, using the identities [41] 



/•oo 
^0 



{t + a-iy. (f 



j+a 



a-1 

r=0 



(0" 



a - 1\ 

r ){t + r)\ 



for integers a and t, and [36] 



/•oo 

/ a;*e-''^da; = r(i+l)a-(*+i) 
Jo 



for integer t > 0. 



(80) 



(81) 



□ 



Appendix III 
Proof of Theorem 2 

Proof: We start by noting that at high SNR, (7) approaches 



X;, =log2det (^^Wfe 

where is an m x to complex Wishart matrix given by 

^ ^ , UkRl for Nr < Nt 
^ H[,Hfc for Nr> Nt 



(82) 



(83) 



Substituting (82) into (47) and using (57), we write the 
variance of the MIMO-OFDM mutual information at high 
SNR as follows 

> N-l 



d=i 



X E 
1 



log2 det ( ;^Wo ) logj det 



Nt 



-E 



log2det{ ^Wo 



E' 



log2det( ^Wo 



(84) 



Noting that 



log2det(Wo) (85) 



log.det(^^Woj=TOlog2(^^^ 

we apply some simple algebra to (84) and find that the terms 
involving 7 cancel perfectly, leaving 



^[log2 det(Wo)log2 det (W^ 



1 



-E 



(log2det(Wo))^ 



N 

i?2[log2det (Wo) 



(86) 



Since Wq is a complex Wishart matrix, we invoke results from 
[3] to give 

m— 1 

E [log2 det (Wo)] = log2(e) ^ i>{n - t) (87) 



t=o 



E 



(log2det(Wo))' ) 



= (log2(e))M EV''(n-i)+E^("-^) ■ 

We now consider the remaining expectation 
E [log2 det (Wo) log2 det (W<j)] in (86). For the exti-eme 
cases of pd ~ Q and pd = 1, this is directly obtained from 
(87) and (88) respectively. The main challenge is to obtain a 
closed-form finite sum expression for < \pd\ < 1. 

We start by following the same procedure as used in (50)- 
(60) in the proof of Theorem 1, which yields 



£;[log2det (Wo)log2det (W^)] 



\Pd\ 



EEdet(D,,,) (89) 



r™(n)r^(TO) (1 - |pd|2) 

for < \pd\ < 1, where 5^,^ is an to x m matrix with entries 
corresponding to (60), but with the «(•) functions replaced 
with 



a{x) = log2(a;) 



(90) 



We now evaluate the integrals for the elements of T>r,s 
corresponding to b{i,j) and c{i,j) in (60), using the identity 
[36, eq. (4.352.1)] 



69 



(V^g) -ln(6)) . g > 0, 6 > (91) 



This gives 



log2(e)r(T + j)|Prfr 



E 

t=0 



J- A ( \pd? VT{u)i^{u) 



t J {r + ty. 



(92) 



and 



c(i,j) 



log2(e)r(T + *)|prfr 



E 



(i-iPdP)- 

i-A f \pd\^ Yt{v)4'{v) 



t J \i-\pd\y (T + t) 



(93) 



To evaluate the remaining integrals in D^.s, i.e. for the 
elements a{i,j), we use (64) and (91) to obtain 



a(z,j) = (log2(e))2|pdr(l-|prf| 



(94) 



E 



\Pd\''r{u)T{v) 

tl (r + t)! 



{H{u-l) + h{pd)) 



X {H{v - 1) + h{pd)) . (95) 
Next we use (92)-(94) in (89), and perform some basic 
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3d I ) V^OO 



oo \pa\^^r{u}r{v){h{pd)+H(u-l)){h(pa)+H{v-l)) 
t!(T+t)! 



for i =/= r, j =/= s 

for i = r, j =/= s 

for i =/= r, j = s 

for i = r, j = s 



(97) 



simplifications to obtain 

i;[log2det (Wo)log2det (W<j)] 

n / \\2 mm 

(log2(e))^ 



^^det(a,,(prf)) (96) 



r™(n)r„(m) 

where Cr,s{Pd) is an m x m matrix with element 
given by (97) at the top of the next page. The expression 
(24) follows by using (96), (88), and (87) in (86). 

To complete the proof we must express the infinite summa- 
tion in (97) in the simplified finite-sum form of (27). This sim- 
plification requires significant algebraic manipulations, which 
we now detail. Start by recalling the definitions u = t + r + i 
and V = t + T + j, and writing the infinite sum in (97) as 
follows 

[c^APd)).= ^]7ljf2 sM) (98) 



where 



, A^x'{T + t + i-iy.{T + t + i-iy. 
^ ^ WTt)i 

t=0 ^ ' 

X {h{V^) + H{T + t + t-l)) 

X {h{V^) + H{T + t + j -1)) . (99) 

Note that the series (99), and those that follow below, are 
convergent for |a;| < 1 (a condition which holds in (98)). 
Now, (99) can be written as 

S{x) = h'^{y^)Si{l,l,x) + h{y^) (^Si{H{i),l,x) 

+ Si{l,H{j),x)^ +Si{H{i),H{j),x) (100) 

where 

Sl{h{l)j2{j),x) 

A x*(t + t + t- 1)!(t + t + j-iy. 



E 



ti{T + ty. 



X /i(r + < + z-l)/2(r + i + j-l) (101) 

for arbitrary functions /i and /2. We now consider each of 
the infinite sums in (100) in turn. 

First consider Si{l,l,x). Following a similar general ap- 
proach to that used in [42], we perform the following sequence 
of operations'* 



Si{l,l,x) = J2 



t=o 



x\T + t + i-iy.{T + t + j - 1)! 

t\{T + ty 



da;^+*- 


- 1 




1 




~i 




1 


da;'^+*- 


-1 




1 


dx^+*" 


-1 




1 


dx^+«- 


-1 




1 


dx^+*- 


-1 



oo 



{r + ty 



OO 



x*+'-^{t + 3 - 1)! 



t=o 



dxJ- 



(102) 



Via application of the Leibnitz formula, it can be shown that 

Si{lXx) 



E 



T + i-i\{i~iy d"+-'+''-i ^ ; 

^ • T+j + b-l Z^^ 



T + b J bl dx'^+j+'' 



b=o ^ ' \t=o y 

h=0 ^ ^ ^ ^ \t=0 



(103) 



Now noting that 



t=o 

with derivatives 



oo 



d"^ 



dx'^ \ / (1 - 



E^ 



|x| < 1 (104) 



(105) 



we can write (103) as follows 

fi - xY+i ^ 

6=0 

r(r + i)x^-i — 



r(T + i) 4:^A-i\/ \'' (r + j + fc- 1)! 



(l-a;)^+J^V ^ /\l-a;y (r + &)! 



(1 — x)^ \ b I \ X 

(l-a;)^ 



(106) 



'^Note that for this particular case, a finite expression could be also found by 
directly matching the infinite series to a hypergeometric function, and using 
associated identities; something which cannot be done in the other cases. S\(H{%) ^ 1, x) 



Now consider Si{H{i), l,x). Following the same sequence 
of operations as in (102) and (103), we find that 
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, b J (t + by. dx^+'+'> , 



(107) 



Now we use [43] 



Yx'Hit)^ ^''^^ |x|<l (108) 

1 — a; 

and the corresponding derivatives 



(109) 



to write (107) as follows 



5i(i7(z),l,x) 



^(t + j) fj -i\ f X y{T + i + b-i) 



(1 - xY+' ^ 



X (i7(T + i + 6 - 1) - ln(l - x)) 



c^-ir(r + j);(-4/j-l\ /I 



E 



X (i?(z-l-6)-ln(l-x)) 



(r + J - 



= _ iVi,j{H{z- -ln{l- x)r],j{l,^/^)) . 

(110) 

Now consider Si{l, H{j),x). Using exactly the same ap- 
proach as for Si{H{i), 1, x), we obtain 

5i(l,i?(j),x) 

= (VjA^iz- -ln(l -x)?7j-,(l,Va:)) ■ 

(111) 

Finally consider Si{H{i), H{j),x). We follow the same 
sequence of operations as in (102) and (103). In this case 
it is convenient to take the successive derivatives based on 
the order of i and j. In particular, with i' — max(i,j) and 
j' — min(i, j), we obtain 



S,iH{i),H{j),x) 



ir+i'-l 



X 



'-1 



i'-l 



6=0 



i' - 1\ x'' d^+j'+b-i 



(112) 



where 



t=0 

In this case, in contrast to the previous summations in (104) 
and (108), the infinite summation in (113) cannot be directly 



expressed in a finite form. To evaluate this series in finite form, 
we start by using (29) to write 



q=l 



t=l 
00 



Y,^'H{tr+ss{x) 

t=i 

Li,(l -x) +ln2(l -s) 
1-x 



53 (x) (114) 



where Lj, (•) is the dilogarithm function [37, eq. (27.7.1)], and 
Sz{-) is given by 



'53(x) ^ E E 



'-^ ^ x'H{t) 



t + q 



(115) 



q=l t=l 

Note that the last line in (114) followed by using an identity 
from [43]. We now manipulate S^{-) as follows 



'53(.) = E-E^^^ 

1 2:9 f-f t + q 



(116) 



Ei,E/-*^''-^^Wd- 

9=1 t=l"' 



13=1 " \t=i 



E-:^/"^'V"^^r^^d.. (117) 



9=1 

For g > 1, consider 



xi 



l-x 



x'i-^ 
l-x 



so therefore 



l-x 

9-1 



l-x ^ 



>1, (118) 



■u=l 



'53(-)-E^E/-""mi 



x)dx 



q=2 v=l 



E:^/'^d.. (119) 



9=1 



a;'^ / 1 — a: 



Using [36, Eq. 2.729]^ 



j 2/"nn(l-y)dy=^(^(y'"+i 



-l)ln(l-y) 



E] ^ ^ — ^ I + const 



fc=i 

and noting that 

ln(l - x) 



m- k + 2 



(120) 



dx = — ln(l — x) — ln(l — a:;)da; 
1 — a; / dx 



'There is a missing (—1) factor in this reference. 
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ln2(l - x) 



const 



(121) 



we can now express 53(2;) in finite form as follows 



q=2 v=l \ 



y - 



2 ^ a;9 

9=1 



(122) 



Note that it can be easily verified, using (116), that the 
integration constant generated in going from (119) to (122) 
is zero. After much algebraic manipulation, it can be shown 
that (122) reduces to 

ln(l -x) ln(l - x) 

q=l 



i'-j'-l 

E 

q=l 



ln(l - x)H{i' - f ~ q) 



ln(l - x)H{q) 

H{r + q-\)-H{r~\) 

xi ' 



, 4 -J - 

- y 



(123) 



Now substituting (123) into (114) we can express 52(2;) as 
the finite sum 



S2{X) 



Li,(l-j;)+ln2(l-a;) 

1 -X 

ln(l — x) ^ — x) 

2 ^ 

9=1 

''^"^ - x)H{i' - f - q) ln(l - x)H{q) 

q=l 



xi 



.9+1 



- E 

j.q ^ 



H{r + q- 1) -H{r- 1) 



(124) 



The corresponding derivatives can be obtained after tedious 
algebra as follows 



dx 



(125) 



where 



e.(r) = Li,(l - x) + ln\l -x)- 2H{r) ln(l - x) 
/2ff(fe-l)-/i,b_i(a 



6=1 ^ 

1 ^ find -gM|) 



9=1 
5-1 



(126) 

where 6 = i' — f , and recall that Li2(-) is the dilogarithm 



function [37, eq. (27.7.1)]. Also, K{-) is a constant given by 
A-M^i: ^" + '-''-^"-" . (.27, 



and 



fqA^)^y^^^-^XqA^)Hl-x) (128) 



where 



Substituting (125) into (112) we obtain 

Si{H{i),H{j),x) 

r(T + i') fi' -i\ f X 



(1 - xY+i' ^ \ b J \ l-x 

^ ' b=Q 



d-xy+i' ^ 



X (t + z'-6)j-'-i^:.(z-5 - 1) 



(1 -x)^ 



(130) 



Finally, substituting (130), (1 1 1), (1 10) and (106) into (100), 
and then combining with (98) and simplifying, we obtain the 
desired finite-sum expression in (27). 

□ 



Appendix IV 
Proof of Theorem 3 

Proof: We start by following [44, 45] and applying a first- 
order Taylor approximation to (7) near 7 = to give 



^/c«log2(e)^tr(HfcHi 



(131) 



Note that, as also mentioned in [44,45], we emphasize that 
this result is only accurate for the low SNR regime; in general, 
requiring that the condition || (7/A^t)HfeHj^.|| < 1 is satisfied. 

Now, substituting (131) into (47) and using (57), we write 
the variance of the MIMO-OFDM mutual information at low 
SNR as follows 

/ \ 2 

VarO(Xofdm) = (log2(e))' ' ^ 



7V2 



y{N~ d)E [tr (HoHJ) tr (HrfH^ 



d=l 



+ [tr^ (Hfl,,HLt)] - [tr {n^.tK.^ 

(132) 

where Hfiat is a flat-fading i.i.d. Rayleigh fading channel 
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matrix. From [7], we have the following results 



E 
E 



tr- 



= NrNt{l + NrNt) 



(133) 
(134) 



For the remaining expectation in (132) we write 



E 



tr ( HoHi ) tr (^H^H;, 



1=1 i=i fe=i 1=1 

= EE^O(Ho)^,,f m,)tj'] + {NrN^y - NrNt 

i=l i = l 

(135) 

where the second line followed by noting that 

i?[|(Ho).,,f|(Hd)fe,,|2] = 1 for all ^ {k,£). 

Now using (5), it can be easily shown that 

£;[|(Ho).^,f|(H,),,,f] 

= [|(H,),j4] + (1 - Ip.Hi? [KEK^n 



Substituting (136) into (135) we find that 



(136) 



E 



tr ( HqHJ ) tr ( HwH 



(137) 



The theorem now follows by substituting (137), (134), and 
(133) into (132) and then performing some basic simplifica- 
tions. 

□ 
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